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Abstract 

Perturbations in Prigogine's cosmological model with particle cre- 
ation are investigated in the framework of Newtonian gravity theory. 

The necessary conditions for the density contrast to reach the non- 
linear regime which is necessary to guarantee structure formation are 
obtained. 

The upper limit of the age of the universe is investigated for the 
perturbations obtained to find them compatible with the observed 
data from COBE(Cosmic Background Explorer satellite). 

We compare the particle horizon and the Hubble sphere for to 
guarantee that the perturbation is in observable universe. 



1 Introduction 



This paper reviews some aspect of our knowledge of the gravitational the- 
ory of fluctuations of density in homogeneous and isotropic models of the 
universe. The evolution of density perturbation is evaluated for Prigogine's 
particle creation model for irreversible process with the framework of Newto- 
nian gravity theory. In this framework various works ^ - [0 study the time 
dependence for the density contrast. We will find some simple conditions 
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on the rate of expansion which permit hydrodynamic perturbation to grow 
rapidally with time. 

Here, we just assume for compatibihty with the cosmological principle, 
that any fluctuation present must have an amplitude which decreases with 
length scales or equivalently, mass scale. 

The importance of the reach of the non linear regime for the density con- 
trast is relacioned with the age of the structures in the universe. Considering 
that the structures in the universe to form "from top down" , we can use the 
period when 5+ ~ 1 and the age of the universe to estimate the age for the 
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2 Prigogine's Model 

The universe has a considerable entropy content, mainly in the form of black- 
body radiation. Therefore, Einstein equations are purely adiabatic and re- 
versible implying in the raising of the question: What is the origin of cosmo- 
logical entropy? 

Prigogine proposed a new type of cosmological history which includes 
large scale entropy production. These cosmologies are based on the reinter- 
pretation of the matter-energy tensor in Einstein's equations and applied to 
the homogeneous and isotropic universe, namely 

ds'^ = dt^ - R{tfda'^ , (1) 
R=R(t) only. The Einstein's field equations 

Gfj,u = kT^j^y (2) 

are used together with the adiabatic transformation in any system for the 
perfect fluid, namely 

T'"' = {p + P)u''u'' + Pg^"' , (3) 

where 

P = P, + Pth. (4) 

Pth is the thermodynamic pressure and Pc corresponds to the matter creation 
that is defined by |jl3|, namely: 

P^ = -P±3ll^n + ^Hn). (5) 
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The process of particle creation is to be consider adiabatic and pressure 
creation is a kind of viscosity pressure that is present in the energy momen- 
tum tensor representing the gravitational fluid. According to the second law 
of thermodynamic the particle number variations admitted are such that 

dN - d{nV) > . (6) 

This implies that, in the presence of matter creation, the usual Einstein's 
equation for (1) is: 

= 3//^ + (7) 



p^-m{p + Pth), (8) 

to become 

kp^3H' + ^, (9) 

and 

71 

P^-{p + PtH). (10) 

n 

In order to exemplify the theory mentioned above, they consider the ther- 
modynamical pressure as zero and the source of the matter creation propor- 
tional to H'^. Thus 

ri + ne = * (11) 

to become 

(12) 

with a > 0. This leads to 

i?(t) = {l + c(e^-l)}3, (13) 

where 

_9_M^i (14) 
kMa^ 3 ' ^ ' 

M is the mass of particle produced and no is the initial number density of 

particles. 
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The evolution of the scale factor is similar to the models with cosmological 
term [14] and have the pattern demonstrate in fig.l. 




Note that the recent measurement of distant supernova when they are 



used as standard candles |jT5[ [[l^] indicates an accelerating universe, conse- 
quently the desaceleration parameter will establish restrictions on the cre- 
ation parameters of Prigogine's model. Namely 



RR 



-{-2c + 3(c- l)e( — —^^}. 
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(15) 



For an accelerating universe q must be negative. This condition imposes 
restrictions for the time coordinate which is: 



t > 
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akM 3c-l 



(16) 



If we wish to have a positive time , consequently c > 3. Therefore, if we don't 
truncate, arbitrarily second Abramo and Lima [1^, the time coordinate at 
t = 0, then the only restriction is given by the condition (16). 



3 Differential equation for contrast evolution 

The fundamental hydrodynamical equations that describe the motion of the 
cosmic fluid here considered are 

i^)r + Vripu) = <il (18) 

V^$ = 47rG^. (19) 

The sub-script r refers to the proper distance referential. 
It should be observed that the particle creation has modified Euler's equa- 
tion [TH], namely 



VP, = ^{up-u) . (20) 

It will be assumed that the velocity of the created particles Up is the 
same as the already existing ones. Equations (17)- (19) are respectively, the 
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momentum conservation law (Euler's equation), the continuity equation and 
Poisson's modified equation. 

Here u is the velocity of a fiuid element, p is the fiuid mass density, $ is 
the gravitational potential ,P is the total pressure . The fiuid we consider 
here is a cold one, such that Pth = and Pc = Pdf) only. 

We now introduce the comoving coordinate x related to the proper coor- 
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f=R{t)x, (21) 

where R{t) is the expansion factor. Equation (21) is a change of variables 
from proper locally Minkowskian coordinates (r, t) to expanding coordinates 
{x, t) comoving in the background model. In these latter coordinates the 
proper velocity of a particle relative to the origin is 

u = ax + xa = ax + v{x, t). (22) 

If a small perturbation is placed on the background, first order corrections 
appear in the velocity, density and potential. We will call these v{x, t), 6{x, t) 
and 0(x, t). We shall assume that these quantities are small such that 5 «1 
and V « u. Peebles |jT^ shows that the potential in the new coordinates 
become 

$ = - ]^RRx'^ . (23) 

We shall not detain to derive equations (17) - (19) in the new coordinates. 
We list them only and refer the reader to [p!9| 



V^0 = 47rG'i?2(p + 3P) (24) 



R_.U..Vu+ — ^^. (26) 

Expanding p , u and (j) pertubatively and keeping only the first order terms 
in the equations above gives the linearized equations. The perturbations are 
defined as follows 
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p = Po{t){l + 6{x,t)) 



(27) 



u^Rx + v (28) 
$ = - ^RRx^ = + '^ttGpqR^x'^ + 2PcR^x\G , (29) 

where 

\v\,\(b\,\6\«l. (30) 
The field equation of zero order has been used, namely 

— = -47rGpo-127rG'Pe- (31) 
The linearized perturbative equations become 

= "^ttGR^poS (32) 

+ = 

V.i; = -it:[^ + *-]. (34) 
dt po 

From these equations we can obtain a second order difi:erential equation 
for the linearized density contrast, S, by taking divergence of equation (33) 
together with (32) and (34) we get 

■■ R ■ R"^ d 

5 + {2^ + -}S - {AttGpo - 2^- - ^(^/po)}S = . (35) 
R Po R Po at 

In order to integrate equation (35), it will be convenient to change variable 
from t to U such that 

R{t) = m . (36) 
After some algebra equation (35) is transformed into 

U\U - a)S" + {^{U-a)U + {h + l)U^}S'-{'^{U-a)- ^hU}S = , (37) 
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where ^3 = ]g and a = 1 — c. The derivative is taken with respect to 
the variable U . By integrating the above equation we obtain the following 
hypergeometric functions 

5^ = U^r\Y{a,,0^.r,'iM/a)] (38) 
52 = U-\Y{a2,P2,l2,U/a)}, (39) 

where 



8 + 3A;3 + ^J^kl 


- 24A;3 + 40 
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8 + 3A;3 - ^JQkl 


- 24A;3 + 40 
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71 = 8/3 




72 = -2/3 


+3A;3 - 2 + y/gfci 


- 24A;3 + 40 
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+3A;3 - 2 - yj%kl - 


- 24A;3 + 40 
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(40) 



A = (41) 

(42) 
(43) 



«2 = (44) 

+3A;3 - 2 - a/qH - 24A;3 + 40 
/32 = ^ . (45) 

It follows that 5{U) represents a family of solutions of equation (37), labeled 
by the parameters ai, /3i,q;2 and /32- It might be noted that the behavior of 
5 in this linear approximation depends only on the initial values of 5 and 



4 Growing and decaying modes 

In order to have an insight about the nature of the solution we consider the 
case where A;3 = 1 and c > 1, such that the hypergeometric solutions and 
5_ become 



5_ = 



Cl 

U{U - a) 
U-a ' 



(46) 
(47) 
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where c\ and C2 are integration constants, 5+ is the growing mode and 5_ is 
the decaying mode. 

Although the decreasing mode can be important in some circumstances, 
we shall hereafter mainly deal only with the increasing mode. It is respons- 
able for the formation of the cosmic structure in the gravitational instability 
picture. Beyond this the universe would not have been homogeneous in the 
past if we consider the coefficient of the decaying mode different of zero. 

Substituting (13) in (17), we can express the mass density contrast as a 
function of the parameter r, namely 



5+ = w\e-^ - (1 - -)e 



-3t1 



(48) 



with w = C2C3 and 



akMt 
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The profile for the 5+ is given below 



(49) 
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Figure 2: Evohition of the growing mode for the density contrast when c — 4. 
For mathematical convenience we take w — 2. 
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and for c = 10 by 



0.2 0.4 0.6 



t 



1.2 1.4 1.6 l.i 



Figure 3: Evolution of the growing mode for the density contrast when c 
10. For mathematical convenience we take w = 2. 



5 The reach of the non hnear regime for the 
density contrast 

In order to evaluate the time for the perturbation to reach the non-linear 
regime for the density contrast in our model, we use the Sachs- Wolf relation 
together with equation (48). Thus 

-7^ oc — 50 

is applicable for the considered models. 



It is known that 20 



(^W< 0(10-2 -10-3). (51) 

If we insert these values in our model, that is, in expression (48) one gets for 
T the following 

Tdec ~ 0(10-^ - 10-5) , (52) 
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where Tdec corresponds to the parameter at the decouphng time. 

It seems reasonable to suppose that the present strong condensations 
grew out of small disturbances so that a necessary (if perhaps not sufficient) 
condition for their formation is that the perturbation, 5+ , calculated in linear 
stability should have become of order unity at some time before the present 

Consequently 5+ ~ 1 implies that 



trl ~ 10-' (53) 



in our model. 

We construct the quocient 

Q = ^ = ^ (54) 

"^dec tdec 

for to determinate the time for the perturbations reach the non linear regime 
(^rl)- Using (53) and (54) we obtain 

g = ^ ^ 0(10^ - 10^) . (55) 

'Tdec 

If we insert these values into (54) and take in account that the decoupling 
occurred at t ~ lO^ys one finds 

0(10^-10'°). (56) 



6 The horizon problem 

Another crucial problem is to verify the horizon problem for Prigogine's 
model. To this end, one should point out the distinction between the cos- 
mological particle horizon and the Hubble sphere, or speed of light sphere, 
Rci which is simply defined to be the distance from O of an object moving 
with the cosmological expansion and the velocity of light with respect to O 
. This can be seen very easily to be 
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The cosmological particle horizon at time t is given by 



Rh = R(t) / Rity^dt . 
Jo 

For our model, this is given by 



(58) 



R^ = - 6- [7 3 In ^2 ^ 2V3arctan— 



+ 2v^arctan 



V3,6i -2, 



f/ is given by (36), b = c — 1 and c is given by (14). Due to the analytical 
difficulties of comparing Rc and Rh we plot their graphs (fig.4). The horizon 
problem is solved when c is large (that is when the initial number of particles 
per volume is large). 



c=2 



c=4 
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c=100 



c=1000 



c=1000000 




Figure 4: Cosmological particle horizon(red) and Hubble sphere (green) vs. 
T for various values of c. The particle horizon intercept the Hubble sphere 
for: c = 2 ^ r ~ 0.38; c = 4 ^ r ~ 0.27; c = 10 ^ r ~ 0.16; c = 100 ^ 
r ~ 3 * 10-2; c = 1000 ^ r ~ 3 * 10"^; c = 1000000 ^ r ~ 3 * 10"^ . 



7 The Age of the universe 

We now turn our attention to evaluate the characteristic time scale for the 
age of the universe with the ultimate aim of determining tp, the time that 
elapsed from the beginning of the universe until now. 

The recent results from the Cosmic Background Explorer (COBE) strongly 
supports the gravitational instability theory of structure formation and can 
be interpreted as providing evidence for the existence of small inhomogeneities 
generated in the early era. For the Priogogine's model the perturbation is 
given by 

R{t) 
^+ = ^ 

where 



(60) 



(61) 



With the aid of equation (60) it can be written as 



fjT^, = (1 + ^.ec)e-'^(*-*-) , (62) 
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where Zdec is the decouphng redshift. 

By assuming that during the matter dominated era particle production 
are restricted to non relativistic ones and its amphtude is fixed at decouphng 
time (COBE roughly fix this amplitude), we have 

and 

S+{tdec, P) = 5+iUec, P = 0) , (64) 

where (Ts{p) is the root-mean-square mass fiuctuation in sphere of radius 
8h~^MPc. The recent COBE measurements indicate that (j8(/3 = 0) ~ 1 for 
the standard cold dark matter [?]. 

By requiring that structures form not to late, it is safe to assume asiP) > 
I Using this condition into equation (63) implies that 

e-/3{tp-w) > 1 _ (65) 

If we assume that tp » t^ec, a superior limit for the age of the universe 
may be presented as 

6.6 

okM 

which becomes in term of the parameter r, using (49), as 

Tp < 0.55 . (67) 

In order to give a better idea on the superior limit of the age of the 
universe, we follow the same procedure used in the last section. We construct 
the quotient 

Q = ^. (68) 

Tdec 

Fixing Tp ~ 0.55 and Tdec ~ 10~^ — 10~^, one gets 

g -0(10^-10^). (69) 

Once again if we take for the decoupling time tp ~ 10^ ys one has then 

tp {10''- 10'') ys. (70) 

This result indicates that for the Prigogine's model, the age of the universe 
is superior to the age of the structure formation and have an interval which 
agrees with the actual estimation. 



tp < -r-^ , (66) 
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8 Conclusions 

In summary, in this paper we analyze the growth of density perturbation in 
Newtonian cosmological model for Prigogine's type model to find a class of 
hypergeometric solutions. 

A special case is studied to obtain an increasing and a decreasing mode 
analogous to Friedman model. 

An estimate of the time interval for the perturbation to reach the non- 
linear regime is obtained. 

An upper limit of the age of the universe is obtained by making use of 
Cobes' data. 

The horizon problem was carefully explained to show that for big values 
of the parameter c the space-time is connected casually. 
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